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Abstract 

In a finite temperature Thomas-Fermi framework, we calculate density distribu- 
tions of hot nuclei enclosed in a freeze-out volume of few times the normal nuclear 
volume and then construct the caloric curve, with and without inclusion of ra- 
dial collective flow. In both cases, the calculated specific heats Cy show a peaked 
structure signalling a liquid-gas phase transition. Without flow, the caloric curve 
indicates a continuous phase transition whereas with inclusion of flow, the transition 
is very sharp. In the latter case, the nucleus undergoes a shape change to a bubble 
from a diffuse sphere at the transition temperature. 

1. INTRODUCTION 

The equation of state (EOS) of nuclear matter constructed with realistic effective in- 
teractions shows that the nuclear liquid coexists in phase-equilibrium with the vapour 
surrounding it below a critical temperature ~ 15 — 20 MeV. Searching for signals for 
liquid-gas phase transition in nuclear matter or in finite nuclei has become one of the 
fundamental issues in medium and higher energy nuclear collisions in recent times after 
the initial experimental observations of a power law behaviour in the mass or charge dis- 
tribution of nuclear fragments in proton and heavy ion induced reactions which is usually 
taken to be indicative of a critical behaviour 0. On the other hand, different models 
of liquid-gas phase equilibrium in nuclei 0, ^ show that there is a limiting temperature 
beyond which the nucleus is unstable and the nuclear liquid and the surrounding gas can 
not coexist in thermodynamic equilibrium. The limitation in temperature also seems to 
be an experimentally supported observation [Q . A possible connection between the phase 
transition temperature and the limiting temperature is therefore called for. Normally 
phase transitions are signalled by peaks in the specific heat C„ at constant volume as 
temperature increases. The microcanonical algorithm of Gross designed for multifrag- 
mentation calculations show such a peaked structure for C^,, so also the calculations by 
the Copenhagen school in a canonical description. Recent calculations by Dasgupta 
et al. [0] in the lattice gas model for fragmentation also show such a peak for the specific 
heat. These calculations however do not allow for any detailed study of the changes in 
the density distribution at criticality. 



Renewed interest in the liquid-gas phase transition has been fueled further by the re- 
cent experimental observation P] that in Au+Au collisions at 600A MeV, the temperature 
in the caloric curve remains almost constant in the excitation energy range of ~ 4 — 10 
MeV per nucleon. This is reminiscent of a first order phase transition and the width in the 
excitation energy at constant temperature may be called the latent heat for liquid- vapour 
phase transition. Another set of experimental data with Au on C is reported recently 
[P] at a little higher energy of lA GeV where the caloric curve looks somewhat different; 
here around 6-7 MeV the temperature rises slowly but steadily with excitation energy 
indicative of a continuous phase transition. These two different behaviours in the caloric 
curve possibly reflect subtle changes in the physical process as the bombarding energy 
increases or as the colliding mass changes. In this paper, we propose an explanation for 
them in a finite temperature Thomas- Fermi theory (FTTF). Part of the work has been 
pubhshed recently . 



2. THEORETICAL FRAMEWORK 

In the following, we give a brief outline for determining the density profiles of hot 
nuclei in the FTTF approximation starting from an effective nucleon-nucleon interaction. 

2.1. The interaction 

The energy density of the nucleus is calculated with a Seyler-Blanchard type mo- 
mentum and density dependent finite range two-body effective interaction It is given 
by 

'Veff{r,p,p) = C/,«K(r,p) +V2{r,p)] (1) 
Vi = -(1 -pV^^)/(n,^2), 
V2 = d'[p,in)+p2ir2)rfin,f2) (2) 

with 

p-\fl-f2\ jfj 

m,r,)= I . (3) 

I ri - r2 I /a 

Here a is the range of the interaction, h denotes its strength of repulsion in the momentum 
dependence, r = |ri — and p = |pi — P2I are the relative distance and momenta of the 
two interacting nucleons. The subscripts / and u in the strength C refer to like pair (n-n 
or p-p) or unlike pair (n-p) interaction respectively, d and n are measures of the strength 
of the density dependence of the interaction and pi and p2 are the densities at the sites 



of the two nucleons. The values of the potential parameters are given in ||Tl 
2.2. Self-consistent density profile and the caloric curve 



by 



The interaction energy density for a finite nucleus with the interaction chosen is given 
2 If 



{Cin^{r',p!) + Cun_r{r' ,p)}n^{f,p)\ dfdpd^, (4) 
and the kinetic energy density is given by 

^('^) = ^Y.J ^n^ir,p)dp. (5) 



Here r is the isospin index, nr{r,p) the occupation probabihty and the nucleon mass. 
The coulomb interaction energy density is taken to be the sum of the direct and exchange 
term. 

The occupation probabihty is obtained by minimising the thermodynamic potential 

G = E -TS -J2firAr, (6) 

T 

where E and S are total energy and entropy of the system at temperature T and A^. the 
number of neutrons or protons, fir being their chemical potential. In the local density 
approximation, the occupation probability can be written as 

2 

nr{r,p) = [1 + ^Mi^^^ + V!\r) + ^^(r) + 5.,,K(r) - fir)/T}]-\ (7) 

where m*{r) is the effective mass dependent on the momentum dependent part of the 
single particle potential p^V^^(r) and V^, V^, V^, Vc etc are potential terms expressions of 



which in terms of density are given in Ref.[ll|. The density is the momentum integral of 



the occupancy. The total energy density is then written as 

e(r) = Pr{r)[Th/2{Vr{r))/Ji,2Mrm - mX{r)V^{r)) + + W^))], (8) 

where Jk are the Fermi integrals and the fugacity riT-{r) is defined as 

^r{r) = W - V^{r) - V^{t) - 5r,V,{r)]/T. (9) 

The total energy per particle E{T) is obtained from the energy density. The excitation 
energy per particle is then calculated as E* = E(T) — E(T = 0). 

The continuum states of a nucleus at a nonzero temperature are occupied with a 
finite probability. The particle density therefore does not vanish at large distances. The 
observables then depend on the size of the box in which the calculations are done. In 
model calculations in heavy ion collisions, it is usually assumed that thermalisation occurs 
in a freeze-out volume significantly larger than the normal nuclear volume. Guided by 
this common practice, we fix a volume, find out the density profiles and obtain the caloric 
curve, the excitation energy as a function of temperature. The specific heat at constant 
volume is then calculated as C„ = {dE* / dT)^. 

Starting from the guess density, one can arrive at a self-consistent density iteratively 



and calculate the physical observables. The details are given in ||TT 
2.3. Inclusion of collective radial flow 

A hot nuclear system created in the laboratory from energetic nuclear collisions may 
be initially compressed. In the decompression phase, a collective radial flow may be gener- 
ated in addition to thermal excitation. An expanding system, in a strict thermodynamic 
sense is not in equilibrium. However, if the time scale involved in the expansion is much 
larger compared to the equilibration times in the expanding complex, the thermodynamic 
equilibrium concept may still be meaningful. In a recent paper. Pal et al. |]I2[ suggested 



that the effect of collective radial flow could be simulated through the inclusion of an 
external negative pressure in the total thermodynamic potential at freeze-out volume. A 
positive uniform external pressure gives rise to compression; similarly a negative pres- 
sure gives rise to an inflationary scenario resulting in the outward radial flow of matter. 
The expanding system can hence be assumed to be under the action of a negative ex- 
ternal pressure Pq. Its magnitude is equal to flow pressure Pj{\Pq\ = Pf), the internal 



pressure exerted by the radially outgoing nucleons at the freezeout surface. The total 
thermodynamic potential of the system at freezeout then is given by 



G = E - TS -Y,lJirAr + PqVL, (10) 

r 

where Pq is the constant external pressure assumed negative and Vt the effective volume. 
It is given as 

^ = ^v^i^^ R = {^-<r'>Yl\ (11) 

Here < >^/^ is the rms radius of the matter density distribution and R is the radius 
of the corresponding uniform density distribution. The occupancy obtained from the 
minimisation of the thermodynamic potential then contains a pressure term and reads as 

n.(f,p-) = [1 + exp{(^^i^ + K°(r) + if (r) + 5.,V;(r) - /i. + P,—Rr^) /T}]-\ (12) 

With the inclusion of flow, the whole set of calculations can then be repeated with an 
effective chemical potential 

^^r" = l^r- Po^-^Rr'. (13) 



It has been shown that the flow pressure -P/(= —Pq) is related to the flow energy as 

Pf = D{vf,T)p{r)ej{r), (14) 

where the quantity D{vf,T) depends weakly on the temperature T and the radially di- 
rected flow velocity Vf and is ~ 4.5 for nucleons and e/(r) the flow energy per nucleon at 
any point within the volume. The total flow energy may be expressed as 



Ef = / p{r)ef{r)df= PfV/D{vj,T). (15) 

Here = V is the size of the box or the freeze-out volume in which the calculations are 
done. 



3. RESULTS AND DISCUSSIONS 



To do calculations in the mass range of experimental interest, we choose two systems, 
namely, ^^°Sm and ^^Kr. We report calculations done in a freeze-out volume V = 8Vq 
where Vq is the normal volume of the nucleus at zero temperature. The freeze-out volume 
chosen here is consistent with most multifragmentation calculations ^ and we have found 
that our results are nearly insensitive to volumes much beyond these chosen confinement 
volumes. Initially, we report calculations without any flow effects. In Figure 1, the proton 
density for ^^"Sm is displayed at four temperatures, T= 5 MeV (long-dashed curve), T= 
9 MeV (dotted curve), T= 9.5 MeV (dash curve) and T= 10 MeV (full curve). With 
increasing temperature, the central density is progressively depleted with thickening of 
the long tail spread at the boundary. Beyond T = 9.5 MeV, the density change is rather 
abrupt and at T close to 10 MeV, the whole system looks like a uniform distribution 
of matter inside the volume. This is shown by a representative density distribution at 
T = 10 MeV. The slight bump at the edge of the density is due to the Coulomb force. 
The caloric curve for the system ^^°Sm is shown in Figure 2 (dotted line). Initially the 
temperature rises faster with excitation energy, then the rise is slower. A kink is observed 
at T RilO MeV, beyond which the excitation energy rises linearly with temperature. The 



0.06 - — — 



150 



Sm 



0.04 



0.02 



0.00 



\ ■ 

\\ \ 

V, 



10 



15 



r (fm) 



Figure 1: The proton density profile for the system ^^^Sm calculated at four temperatures 
as mentioned in the text. 
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Figure 2: The caloric curve for the system ^^^Sm at three pressures. 
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Figure 3: The specific heat per particle plotted as a function of temperature for ^^°Sm. 
The three curves refer to Po = 0.0 (dotted), Pq = -0.05 (dashed) and Pq = -0.1 (full 
curve) respectively. 

slope is then 3/2 which is a reflection of the fact that the nucleons then behave almost 
like a classical gas of particles. In Figure 3, the specific heat at constant volume Cy is 
displayed (dotted fine). It shows a sharp peak at a temperature T ~ 10 MeV. We beheve 
this to be possibly a signal of the liquid-gas phase transition. The abrupt transition in 
the density at this temperature and the sharp fall of to the classical value of 3/2 lend 
credence to our surmise. Collisions between heavy nuclei at high energies may generate a 
modest amount of compression even for far-central impacts. We have therefore repeated 
the calculations taking into account the effect of fiow energy with different values of Pq. If 
Po is given, an estimate of the fiow energy per nucleon can be easily made with the help of 
Equations (14) and (15). For example, if Pq = —0.1 MeV fm~^, the average flow energy 
per nucleon is ~ 1.3 MeV. The caloric curves for Pq = —0.05 (dashed line) and —0.1 
MeV fm~^ (full line) are displayed also in Figure 2. We find that with increase in fiow 
energy, the rise in temperature is slower and when the pressure Pq = —0.1 MeV fm~^, the 
caloric curve shows a plateau at T ~ 5 MeV in the excitation energy range of 5-10 MeV. In 
Figure 3, the corresponding specific heats are displayed. The broad-based peak for no fiow 
goes over to an extremely sharp peak from T ~ 10 MeV to T ~ 5 MeV with increasing 
fiow. Looking at Figures 2 and 3, it is evident that the system signals a liquid-gas phase 
transition at the peak temperatures and that the system moves from a continuous phase 
transition to a sharp first order phase transition. With increase in temperature, it is 
natural to expect that the rms radius of the nucleus in question would increase. That is 
displayed in Figure 4. Interestingly, we find that for cases accompanied with fiow energy, 
there is an extremely sharp increase in the rms radius at the phase transition temperature. 
From Figure 5, we notice that near the phase transition temperature, within an interval 
of AT « 0.1 MeV, the density disribution now undergoes an exotic shape transition to a 
bubble shape at T= 5.3 MeV. This may be attributed to the loss of surface tension at the 
transition temperature; when there is flow, the flow pressure pushes the matter outwards 
which takes a bubble shape because of the constraint of a spherical freeze-out volume. 
The calculations have been repeated for ^^Kr with freeze-out volume 8 times its normal 
volume and all the features reported for ^^°Sm are reproduced except for the fact that the 
transition temperature is shifted up by ^ 0.5 MeV and the peak for the specific heat is 
broader. 
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Figure 4: The proton rms radius as a function of temperature for ^^'^Sm at three pressures. 




Figure 5: The proton density profile for ^™5'm with Pq = -0.1 MeV fm ^ for T = 5.2 
MeV (dashed hne) and T = 5.3 (full hne) MeV. 



4. CONCLUSIONS 



We have calculated the caloric curve and the specific heat of finite nuclei in a self- 
consistent finite temperature Thomas-Fermi theory with a realistic effective nucleon- 
nucleon interaction. The effects of collective radial ffow, if any, are built in the theoretical 
framework. /,From the caloric curve and the specific heat, signals of liquid-gas phase tran- 
sition are obtained for cases both with and without collective fiow. In the case without 
fiow, there appears to be a continuous phase transition at a temperature T ~ 9.5 MeV 

whereas with inclusion of a little flow energy, the phase transition seems to be of flrst 
order at a temperature T ~ 5.0 MeV. The nucleus also then changes its shape from a 
diffuse sphere to a bubble. A qualitative reference to the two sets of experimental data 
referred to in the introduction may not be out of place. In the Au induced reaction on C 
at 1 AGeV, for near peripheral coUisions, there is possibly little or no compression because 
of the small target size. The caloric curve then alludes to a continuous phase transition 
as theoretically obtained and experimentally reported. For the Au on Au reaction at 
600 AMeV, a modest compression and hence radial flow may not be ignored and then a 
sharper phase transition may result as indicated in the experiments. 
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